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Introduction 


Joseph Fourier, whose name forms part of the title of this unit, has already been 
mentioned in this course, in Unit 12, on heat transfer. The present unit is about a 
different part of his work —one which led to results which have become an 
indispensable part of the theoretical tool kit of applied mathematicians, physicists 
and engineers. 


The fundamental idea is that many functions can be represented by the sum of a 
series of sinusoidal terms. In other words, sinusoidal functions may be considered 
to be the building blocks from which other functions may be built up. The point 
of such a procedure is that in many cases a series of sinusoidal terms is easier to 
deal with than the original function. 


The study of the use of sinusoidal functions to represent other functions is called 
Fourier analysis and is the subject of this unit. 


Study guide 


The sections of this unit are intended to be read in the order in which they appear. 
The first four sections all teach new and assessable material. Section 5 does not 
teach any new material but consists of problems on the material of the rest of the 
unit. You should try as many of these as you have time for without looking at the 
solutions. Then compare your answers with the given solutions, and also look at 
the solutions of any problems you did not attempt— you may well find them 
instructive. 


There is a television programme associated with this unit. You should preferably 
have worked through the first three sections before you watch it but, even if you 
have not done so, it is worth watching the programme rather than missing it 
deliberately. There is no television section in the main text of this unit but there 
are notes about the programme in Appendix 3: you should read these notes before 
watching the programme. There is no tape associated with this unit. 


1 What Fourier analysis is used for 


The Fourier method is to express a function as a sum of sinusoidal terms. This 
representation is particularly appropriate when the function is periodic. In this 
section we shall look briefly at periodic functions and then, through examples, see 
why it is useful to be able to express a general periodic function as a sum of 
sinusoids, 


1.1 Periodic functions 


A function is said to be periodic if its graph consists of a continually repeated 
pattern, as in Figure 1. We can express this more formally by saying that a 
function f(r) is periodic when 


S(t) =f(t+nT) for all t, forn = +1, +2, +3, 


Definition (1) implies that the graph of f(t) is unaffected by being shifted a 
distance nT to the right or left. 


(1) 


The positive constant T is called a period of the function. Usually T is taken to be 
the size of the smallest repeating unit in the graph of the function, as in Figure 1, 
but note that we could take any integer multiple of this basic repeating unit as the 
period since this, too, will satisfy the definition (1). 


When we want to apply the Fourier method to a general periodic function we 
shall need to be familiar with the periodic behaviour of our building blocks, the 
sine and cosine functions, For example, the function sin (2t/T) has period T 
because 


_ 2n(t+nT)_ 
sin =s 


2nt +See mid 2nt 
T in T 7} = si 


T° 


On) 
SS et 
period 
(a) 
fO4 
a 
= 
period 
(b) 


Figure 1. Two periodic 
functions 
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so we would expect the function sin (2mt/T) to be useful when we approximate 
functions with period T. 


Question 


2 
Does the function sin an where m is an integer, have T as a period? 
2; 
Does cos? 


Answer 
Yes to both questions, because 


sin Ct +nT)_ si 2nmt Pee 
T = sin{—-~ + 2nm 


. 2mmt 
=sin—_, 


T because mn is an integer, 


and 
og aE +nT) = 2nmt 
Zz 7 T 


We shall see later how these functions are used in constructing Fourier 
representations of functions with period T. 


by a similar argument. 


1.2 An application of Fourier analysis 


Fourier’s idea of using a sum of sinusoidal terms to represent a periodic function 
turns out to be very useful in many different fields. One of these is the response of 
mechanical systems to any arbitrary periodic input. Since you have some 
knowledge of such systems from Units 7, 8 and 24, the following examples and 
exercises are drawn from that field. 


Example 1 

Assume that the system we want to study can be modelled by the lumped- 
parameter spring-mass system with linear damping shown in Figure 2. (This is 
similar to the systems you modelled in Unit 8.) The force, F, is periodic but non- 
sinusoidal. We want to know how the position, x, of the particle, measured from 
its equilibrium position, will vary with time. How can the Fourier representation 
of F be used to find x? 


Solution 
The equation of motion is 


mX(t) + rX(t) + kx(t) = F(t). (2) 


It could well be very difficult to solve this equation with F in its original form, as 
a general periodic function. But now suppose that we can express F as a sum of 
several sinusoidal terms, say 

F=Fy+F,+Fy+0°+Fy 


We can use the fact that Equation (2) is a linear second-order differential equation, 
This means that the superposition principle can be applied. Thus if we can find the 
steady-state outputs, X1, X2, ..., %, that would be produced by the component 
inputs, F,, F3, «++, F,, acting separately, then we know that the steady-state 
output, x(t), produced by their sum, F, will just be the sum of these: 


X(t) = X(t) + X2(t) + -** + x(t). 


Now, since the terms F,, F, ..., F, are sinusoidal, we have reduced the problem 
of finding x(t) to that of finding the output, x,(t), produced by a single sinusoidal 
term, F,(t), What we need to know is the frequency response of the system, that is 
the way in which the output produced by a sinusoidal input of given amplitude 
depends on the frequency. But we do know how to find the frequency response, 
using the methods of Unit 8, 


In this unit I shall use the results we obtained in Unit 8, but with some changes in 
notation, (The notation I use here is common in the sort of engineering 
applications where Fourier methods are used.) In Unit 8 we represented a general 


Strictly speaking these two 
functions have smallest period 
T/m but, since any integer 
multiple of T/m is also a 
period, the functions also have 
T as a period. 


Figure 2, Model system. [p is 
the unstretched length of the 
spring. 


Equation (2) can be derived 
by the methods of Unit 8, 


You met the superposition 
principle in Section 2.4 of 
Unit 6, 


6 MST204 31.1 


sinusoidal input of angular frequency w by Ag cos (wt + oo). Here I shall use the 
input Fo sin wt, which corresponds to @y = —x/2. I shall write the corresponding 
steady-state output as 


x = Asin (wt + ¢), 


¢@ is the phase of the output relative to the input, and is often called the phase shift. 


We want to find x for a given input amplitude Fo and different values of the input 
angular frequency, @. We can do this by investigating the behaviour of the relative 
response, A/F, and the phase shift, @, for varying values of the angular frequency, 
«. This can be represented by the two graphs sketched in Figure 3, which together 
show the frequency response of the system. The exact shapes of the graphs will 
depend on the amount of damping. Notice that in this example @ is negative for 
all values of w. 


Graphs like those in Figure 3 enable us to read off yalues of A/Fy and of ¢ for 
any given angular frequency. So we can find the output, x;, for each sinusoid, F,, 
that makes up the original input. We can then add up all the x, to find the total 
response to the periodic input, F. 


In the above example we used the fact that the system was modelled by a linear 
second-order differential equation, so that we could use the superposition 
principle. More complicated systems may lead to higher-order differential 
equations, but so long as the equations are linear, the superposition principle still 
applies. We define a linear system as a system that can be modelled by a linear 
differential equation. 


We can now generalize the method of Example | to find the response of any linear 
system to a periodic input, F. If we can express F as a sum of sinusoidal terms 
then all we need is a knowledge of the frequency response of the system. This may 
be found analytically, as in Example 1, or experimentally, by observing the 
response of the system to sinusoidal input forces of various frequencies. (Such 
experiments afe standard practice in many practical situations.) 


Example 2 


A linear system is subject to a periodic input, f(t), which can be expressed as the 
sum of two sinusoids as follows: 
S(t) = 2 sin St + 15 sin 101 
where ¢ is time. 
(0) Figure 4 shows graphs of y, = 2 sin St and y, = 1,5 sin 10¢. Plot f(t) and 
determine its period. 
ya 


_— yy = Zsin St 


axe yp = 1.5 sin 10 


7 \ 


0.102 04 04 05 0. 


Figure 4 


In Unit 8 we used the phase 
ne A@, to describe the 
ifference between input and 
output. The phase shift is 
Ad. 


AIF 
0 o 
ry 
0 o 
- 2 
-m 


Figure 3. Frequency response 
for the model system in Figure 
2 


The definition of a linear 
second-order differential 
equation in Unit 6 extends 
naturally to higher orders. 
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(ii) Figure 5 shows the frequency response of a linear system. Derive an 
expression for the response of this system to the input f(t). 
4 
output amplitude 


input amplitude: 


i i 
40 coe 
(a) 
eel 10 20 30 40 30 w 
F TST SHE RAGE Titnhuinhe 
Be Hit Hatt ia 
=a HH it aE 
GHEE: 2H: Heli 
if PSUS ed EEE 
— 40 FESS i 
SEE E if HEHE: 
Hah is : 
oor 1, = 
| ies: 
ft is 
=80° i 
E Subtsbsatasis te HG 
ia HH i 
100° if HE SSEEEIEn aiztit Note: the phase shift is 
measured in degrees. 
(b) 
Figure 5. Frequency response of a linear system. 
Solution 
@ yy 
34 y= yy + yo =2sin Se+ 1.5 sin 10 
2 
y; =2sin5¢ 
- : 
ia yo = L.Ssin 10¢ 
14 ih a 
\ / 


o 0.1 0.2 04 04 be 08 0.9 yo Mi 1.29 13°14 t 
x i i 
t \ / : 
a] Sy 
See 
-24 
=a4] 


Figure 6 
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Figure 6 shows both components, y, and y2, of f(t) and their sum, y, which equals 
S(t). The pericd, T, of f(t) is equal to that of the component 2sin 5t; in other 
words 2n/5 = 1.26. 


(ii) To find the response to f(t) we find the output corresponding to each input 
term separately. 

For the output corresponding to the term 2sin St we need to find the relative 
response and phase shift for the frequency of this input, ic. w = 5. 

From the top graph in Figure 5 we read the relative response for w = 5 to be 


about 0.93. So the output amplitude for the input 2sin 5¢ is 0.93 times the input 
amplitude, i.e. 


093 x2=19. 


From the lower graph in Figure 5 we read the phase shift, @, corresponding to 
«@ = 5 to be about —26°, or —0.45 radians; so the component of the output 
contributed by the input term 2sin Sr is approximately 


1.9 sin (St — 0.45). 


(Remember, from Unit 8, that the frequencies of input and output for a linear 
system subject to a sinusoidal input are equal.) 
Similarly, for the output component corresponding to the input term 1.5sin 10¢, 
we read:off the relative response and phase shift for w = 10. The relative response 
is about 0.72, so the output amplitude at @ = 10 is 0.72 x 1.5 = 1.1. The phase 
shift is about —46°, or —0.8 radians. Hence the output component contributed by 
the term 1.5sin 10r is approximately 

1.1 sin (10¢ — 0.8). 


Consequently, by the superposition principle, the total response of the system to 
J (t) is about 


1,9sin (St — 0.45) + 1.1 sin (10¢ — 0.8). 
The point that is made in part (i) of the solution is a particular example of a more 
general principle: 


If you take the sum of a number of sinusoidal functions of the form A, sin wt, 
Azsin 2ot, Ay sin 3ot, +++, A,sin ket, so that all the frequencies are integral 
multiples of the lowest frequency, then the sum is periodic and the period is equal 
to that of the component with the lowest frequency (often called the fundamental 
angular frequency). In other words the function 


S(t) = A, sin@t + Az sin2ot + Ay sin 3mt + «+++ A, sin kor 
is periodic and its period is T = 2n/o. 


Exercise 1 


Suppose the variation with time r (in seconds) of a voltage at a point in an electrical 
network is given by 


V(t) = 0.01 sin 2r + 0.005 sin 4r + 0.002 sin 82. 
Determine 


(0) the fundamental angular frequency; 
(ii) the period of the voltage. 


[Solution on p. 23] 


Exercise 2 
Hexe 7 shows the frequency response of a linear system. The input of this system is given 
y 
S(t) = sin 4t + 2sin 81 + 3 sin 12¢ 
where t is the time in seconds. 
(i) What is the period of f(t)? 
(ii) Derive an expression for the response, x(t), of the system to the input, f(r), 
[Solution on p. 23] 
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Figure 7 


As I have already said, although Figure 2 shows the kind of model that is used in 

the analysis of mechanical vibrations, and Figures 3 and 7 look like the frequency 

response associated with such a model, you must not conclude that the procedures 

in this unit are useful only in the context of mechanical vibrations, On the 

contrary, they are widely used in such different fields as acoustics, meteorology, 

electrical network theory, and the design of measuring intruments; and this is by 

no means a complete list. In addition, Fourier series are mathematically useful in, See Unit 32, Partial differential 
for example, the solution of differential equations. equations, 


So, in order to decide such different matters as the suitability of a given hi-fi 
amplifier for the reproduction of a particular tape recording or the relationship 
between the reading produced by a measuring instrument and the actual yalue of 
the measured quantity, we need the frequency response of the amplifier or the 
measuring instrument and a procedure for expressing the signal (i.e. the sound or 
other quantity of interest) as a sum of sinusoidal functions. The next section is an 
account of such a procedure. 


Summary of Section 1 
1. A function fis periodic if 
S(t) =f(t + nT) for all ¢, for n= +1, +2, +3,..., 
where T is the period. 
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2. The way that a mechanical system responds to sinusoidal inputs of different 
frequencies is called the frequency response of the system. 


3. A linear system is one that can be modelled by a linear diffierential equation. 
A linear second-order differential equation was defined in Unit 6 —this 
definition extends. naturally to higher-order equations. The superposition 
principle applies to any linear differential equation. 


4. If the frequency response of a linear system is known, the response to any 
periodic input, F, can be found provided that F can be expressed as a sum of 
sinusoidal terms, 


FoF, +Fr+-+ Fy 


The output, x, will then be the sum of the outputs, x;, to each sinusoidal input, 
F,, taken separately. 


5. The sum of a number of sinusoidal terms of different frequencies is a non- 
sinusoidal function. If the frequencies of the sinusoids are integral multiples of 
the lowest one, the fundamental angular frequency, «, say, then their sum will 
be periodic with period 27/w, i.e. the period corresponding to the fundamental 
angular frequency. 


2 Fourier series for periodic functions 


2.1 Introduction 


The aim of this section is to set up a systematic procedure for writing a given 
periodic function in terms of the sum of sinusoidal (or ‘harmonic’) components. 


From the last section you will remember that the sum of a number of terms like 
A, sin wt + Az sin 2ot + As sin 3at + «+> + Aysin kot (1) 


(where Ay, Az, As, .-., A, and @ are constants, with w > 0 and k a positive 
integer), will be a periodic function, f, of r, with period 2x/w which equals the 
smallest period of the term A, sin wt with the lowest, or fundamental, angular 
frequency. What we are aiming at is to start with a periodic function f and end up 
with a series like (1), with all the constants taking known values. Before we get 
down to this, I want to make three adjustments to the series (1), all of them aimed 
at making it more general, that is to say, capable of representing more functions. 


One adjustment arises from the recognition that what we have so far done with 
terms like sin wt and sin 2wt can also be done with terms like cos wt and cos 2wt, 
In other words, you can get a periodic function by adding up a number of terms 
like 


B, coswt + B, cos 2mt + By cos3wt + ++: + By cos ket. (2) 


To take account of this fact I shall represent the series which constitutes our initial 
function by the sum of a number of sine and cosine terms, like this: 


A, sin wt + Azsin2mt + Ay sin3mt + ++ + A, sinker 
+ B, coswt + Bz cos 2mt + B; cos 3mt + ++: + By cos ket. 


This can be written more concisely as 

k 

¥ (A, sin neot + B, cos nat). 

nt 
The second change I want to introduce removes another restriction which arises 
from the fact that the mean value of any function like A sinnet or B cos nwt taken 
over any interval of length T is equal to zero. Any sum of such terms also has a 
zero mean value, like the function in Figure 1(a) of Section 1. Now the procedure 
I am about to show you is not restricted to functions with a zero mean. It can 
quite easily cope with other mean values. For example, the periodic function 
sketched in Figure 1(b) of Section | has a positive mean value and our procedure 
is quite capable of dealing with a case like that. It can also cope with a negative 
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mean value. In order to allow for a non-zero mean, all we need to do is to shift 
the graph of our function up or down by the required amount, which simply 
requires an appropriate constant added to the series, so that our series becomes 


k 
S(t)}=M + ¥ (A, sinnot + B,cos net) (3) 
nat 


where M is the new constant. 


The addition of the term M will not, of course, change the period of the right- 
hand side of Equation (3). As far as frequency response calculations are concerned 
(like those we did in Section 1), M is a cosine term with zero frequency and its 
contribution to the total response is worked out accordingly, as in the example 
below. 


Example 1 
Work out the effect on the response of the system described in Exercise 2 of 
Section 1 of adding a constant term to the input function f(r) so that the new 
input function becomes 

S(Onew = 3 + sin 4t + 2sin Br + 3sin 121, 
Solution 


The new constant will add an extra term to the response. It corresponds to a 
cosine term with zero frequency, From the graph in Figure 7 of Section | we see 
that the amplitude ratio is 1 and the phase shift is zero for @ = 0. So the output 
term corresponding to 3 in the input is 3 x cos(0 + 0) = 3 and the total response 
will now be 

X(Dpew = 3 + 1.2sin (4¢ — 0.19) + 4sin (8¢ — 0.77) + 4.5 sin (121 — 2,3), 


There is one more change I need to make to the series (3) so that it can represent 
a wider class of functions, and that is to extend it to an infinite series, i.e. 


SO=M+ y (A, sin not + B, cos nent). (4) 
mt 


The right-hand side of Equation (4) is known as the Fourier series for f(t). Now, 
we know that this series represents a periodic function with period 2z/a. But what 
we are trying to do is to find this series for any given periodic function. Can this 
be done? And, if it can, how do we find the appropriate values of the 

constants M, A, and B,? 


The answer to the first question, for practical purposes, is “Yes’. (We shall see what 
conditions the given function must satisfy later, in Subsection 2.3.) But for now we 
shall assume that our periodic function can be expressed as a Fourier series and 
get on with finding the constants. 


2.2 Finding the constants 


Let us see how to calculate the constants M, A,, By, A2, Bz, ... in Equation (4) 
from our knowledge of the function /(t) itself; this might be given, for example, in 
the form of a graph. We shall need to use the following integrals. 


Some useful integrals 

In all cases, m and n stand for non-negative integers and « = 2n/T. 
en2 1/2 

[sinnordt=[ — cosnordr=0, #0. (sa) 

J-12 Jr2 
1/2 1/2 

i sin net sin mont dt = J cos nwt cos matdt = 0, n#m. (5b) 
=1/2 =1/2 
7/2 

f sin nwt cos mat dt = 0. (5c) 
-T2 
7/2 Ti2 T 

f sin? not dt = f cos? not dt = zy nO. (Sd) 
-1/2 =1/2 


You can check these integrals 
using the Handbook. 
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Now to find M we integrate every term on both sides of Equation (4) with respect 
to t between the limits — 7/2 and T/2. Hence 


T/2 7/2 T/2 T/2 
S(t)hdt= Madt+ f A, sin wt dt + f B, cost dt 
“7/2 -T/2 -T/2 -T/2 
7/2 1/2 
+f Azsin2ord + f B, cos2wtdt + *--. 
~T/2 -T1/2 


On the right-hand side, every term except the first one is zero according to (5a), so 
that 


7/2 
S(t)dt=T™M 
7) 
1 1/2 
therefore M = 7] , J (t)dt. (6) 


In other words, M is the mean value of f(t). 


We can find A, by multiplying every term on both sides of Equation (4) by sinwt 
as follows: 


S(t)sin@t = Msinwt + A, sin? wt + B, coset sinwt + A; sin 2otsin wt 
+ Bz cos 2tsinwt + °°: 


and now we integrate every one of these terms over ¢ between the limits of — 7/2 
to T/2; 


7/2 T2 1/2 
S(Osinwtdt = Msinerdt + J A, sin? wrdr 
-72 -T2 -T2 
7/2 72 
+ B, cosmtsin wrdt + f A; sin 2tsinqrdt 
-72 -T.2 


1/2 

+ f B, cos 2ewrsin wt dt + +++. 
“1/2 

If we compare all the terms on the right-hand side with (5a), (5b), (5c) and (5d) we 

see that they are all zero except the second one. Hence 


T/2 7/2 op 
Fsinorde= f A, sin? wtdt = A, x — 
72 -T/2 2 


2 7! 
Ay == J (t)sinet dt. 


Pisa 7/2 
Exercise 1 
By multiplying every term on both sides of Equation (4) by coswt and then integrating 


each of the resulting terms over time from — : to z show that 


2 pte 


B, == (Ocoswtdt 
Serie 


where T is the period of f(t). 
[Solution on p. 23] 


All the other coefficients are found in the same way: by multiplying all the terms 
by the sinusoidal function which forms part of the same term as the required 
coefficient and then integrating over time from —T/2 to T/2. It follows, for 
example, that 


7/2 T/2 
J (o)sin Zot dt = f Az sin? 2ot dt 
72 


-T/2 


so that 


2 T/2 
we al f(e)sin 2wrde. 
T J-a2 
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Similarly 


2 1/2 
B,= =f S(t)cos 2er dt and so on. 
T Jen 


Hence for any value of n 


2 rm 
n= = [ J (t)sinneot dt (7) 
T J-42 
and 
1/2 
J (t)cos not dt. (8) 
-1/2 


By following the procedure in the box below, we can now find the Fourier series 
corresponding to some simple functions. 


Procedure 2.2 
To find the constants of the Fourier series 


S(t)=M+ ¥ (A,sinner + B,cos ner) (4) 


nat 
for a periodic function f(r) starting with a graph of f(t) against t. 
1, From the graph determine T, the smallest period of the 


function f(t), and hence the fundamental angular frequency « 
by putting w = 2n/T. 


T/2 
2, Write down the equation for Sf (t)dt in terms of 1. 
“7/2 
1 T/2 
3. Find the constant term M = — S(t)dt. (6) 
T J-rp 


(This is the mean value of f(t) over one complete period.) 


4. Find the coefficients A,, A, As, ... of the sine terms by evaluating 
2 7/2 
pale 


(w as determined in Step (1) above), 


J ()sinneot dt (n= 1, 2,3, ...) (7) 
2 


5, Find the coefficients B,, Bz, Bs, ... of the cosine terms by 
evaluating 
1/2 


Bee J (t)cos nest dt (n=1, 2,3, ...) (8) 
r2 


Ta 


(@ as determined in Step (1) above), 


Example 2 


Derive the Fourier series for the saw-tooth function shown in Figure 1. 
OLY 
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Figure 1. Saw-tooth function. Note that this function is not defined at t= 
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Solution 
First, we find the period, which equals 2x because the graph repeats at intervals of 
2n. From this we find the fundamental angular frequency @ to be | from the 
formula w = 27/T. It follows that the Fourier series has the form 
S()=M + Y (A,sinnt + B, cos nt). 
nt 


Now M is the mean value and is given by 


1 T/2 a h i 
M=— t)dt where > = 7. 
T nih (t) 2 
For t between —7 and 7, f(t) =’t/n (except for the end-points, 1 = —m and t 


which do not affect the value of the integral). 


Therefore 
Wye Ne 
-x J f4- [5] 


=n 


(This was to be expected from the symmetry of the graph.) 
Again 
7/2 


A = J (t)sinneot dt a  sinntdt 
I JA fees nJigt 


1" 
- =f tsin nt dt. 
T do» 
Now, from the tables of standard integrals in the Handbook we find that 
fesin aca = — cost + + sinat +. 
a a 


In our case, a = n and for a definite integral we may omit the constant C, so that 


1 t ee 
y= A ~Leosnt + Jsinm | 
T n n 


=< [-neosne +0 + (—n)c0s(—ne) ~ 0] 


—2ncosnn _ —2cosnm 
ee | eae pea 
Now cosnm = (—1)" so If you have not met this 
before, note that 
Ay =a (a lh ma ei ptt feeant 
“keg EL and so on, 


(Remember that n is any positive integer.) 
So 


2 
A=5 ATs A=— Ag = —, and so on. 
Similarly 
2 7/2 1 e 
B= >[ Fle)cosnondt = % { toos ntdt. 
T Jr To dan 
Again, we find from the Handbook that 


fa i 
fecosarae = tinat a sat +C. 


Again, a = n and C is omitted for a definite integral, so that we have 
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Hence the required Fourier series consists only of sine terms since M = B, = 
therefore 


=F 2 "'si 
fO= & a 1)"*! sinnt 


fran _sin2r | sin3t _ sin4r 
™ 2 3 4 


Notice that the higher the frequency of a term, the lower its amplitude, 


Exercise 2 

How would the Fourier series derived in Example 2 above be affected if 

(i) the height of the ‘saw-tooth’ were doubled, i.e, the value of f(r) increased from —2 
to +2 as ¢ increased from (2m + 1)n to (2m + 3)n (where m is any integer)? 


(ii) the height of the saw-tooth remained the same as in Figure | but the period were 
doubled to 4x? 


[Solution on p. 23) 


2.3 What kind of functions? 


What we have done so far is to show that if a function f(t) with period T can be 
expressed in the form 


S(t\)}=M + ¥ (A,sinnat + B,cos nest) (4) 
wel 


then the coefficients M, A,, B, can be calculated from the formulae (6), (7), (8). To 
make full use of this result we need to know under what conditions a given 
periodic function f(t) can be expressed in the above form. 


There is a theorem which gives a precise result of this kind; in order to state it we 
need two new mathematical concepts: continuous functions and piecewise 
continuous functions. 


A function is said to be continuous if its graph is an unbroken curve, that is if 
there are no gaps in it. For example the sine function is continuous, but the ‘saw- 
tooth’ function of Example 2, illustrated in Figure 1, is not, because of the gaps in 
the graph at r=..., —7, 2, 3m,.... 


A function is said to be piecewise continuous if its graph consists of continuous 
pieces; there must be only a finite number of pieces on any finite part of the 
domain, and the two pieces ending at each gap must have well-defined end-points 
even though these points need not themselves be parts of the graph. For example 
the saw-tooth function illustrated in Figure 1 is piecewise continuous; each piece 
has well-defined end-points even though these are not part of the graph. On the 
other hand the function 1/x is not piecewise continuous because the pieces do not 
have well-defined end-points at the gap where x = 0. 


Now we can state the theorem: 


Theorem 2.3 

If a periodic function f has period T and is piecewise continuous, 
and if its derivative is also piecewise continuous, then it can be 
represented by the Fourier series 


S(t}=M+ y (A, sin nwt + B, cos net) (4) 
=1 


where w = 2n/T. The coefficients are given by 


1 7? 

M=— (tat 6 
T Pe (6) 
2 T2 

A== Sf(t)sinneot dt (7) 
T Jor 
2 Te 

B,== Sf(t)cos net dt. (8) 
T Ja 


The formula (4) is valid for all values of t at which there is no gap 
in the graph. 


The graph of I/x is shown in 
the Handbook, 
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Exercise 3 
Derive the Fourier series for the ‘square wave' function shown in Figure 2. 
1 | 4 I, | i 
| eke | | \ 
1 | | ! \ 
1 | | \ | 
i | \ | \ 
| oa \ | i 
| i | | i 
| | ! | | 
| i 1 1 
1 | i \ \ 
i | | \ | 
| | \ \ | 
\ | | \ i Pe 
32 = 12 a2 32 a 
1 ! | | 1 
| | | | \ 
\ | | | | 
\ | 1 | | 
| \ | 1 \ 
| | | 1 \ 
| i | | | 
1 \ i ' ' 
Figure 2. Square wave function. Note that the function is not defined at t= ..., ~3/2, 
=1/2, 1/2, 3/2)... 
[Solution on p. 23} 
Exercise 4 
Derive the Fourier series for the periodic function whose graph is shown in Figure 3. 
Koa 
: ~ 
4-3-2 = 1 2 3 4 y 


Figure 3 
[Solution on p. 24) 


Exercise 5 
For the function shown in Figure 4, show that M = B, = 0 and derive an expression for A,, 


Figure 4 
[Solution on p. 24) 
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Summary of Section 2 


1, The Fourier series for a function f(t) with period T is 


S()}=M+ y (A, sin not + B,cos nest), 


Fest 
where @ = 2n/T. 
2, The constants in the Fourier series are given by See Procedure 2,2, 

1 

=> t)dt 
TI. aif ( 
2 1/2 
A,== SF (t)sin not dt (n= 1, 2, 3, ...) 
T J-12 
1/2 


J (t)cos net dt (n=1, 2,3, ...). 
a 


3. (i) A function is said to be continuous if its graph is an unbroken curve. 

(ii) A function is said to be piecewise continuous if its graph consists of 
continuous pieces and if there are a finite number of such pieces on any finite part 
of the domain. The end-points of the pieces must be well-defined, 


4. Any periodic function which is piecewise continuous and whose derivative is Theorem 2.3, 
also piecewise continuous can be represented by a Fourier series. 


3 Some outstanding points 


3.1 Odd and even functions 


The object of this subsection is to show you how to use the symmetry of the graph 
of f(t) to save yourself work in calculating the Fourier series. 


In Section 2 we derived the Fourier series (or at least the values of the coefficients) 
for the four functions which are sketched in Figure 1. 


fot (0) 


ff 1 


(a) (b) 


Figure 1. Periodic functions with symmetry about the origin. 


All these four functions have some symmetry about the origin. But the symmetry 
takes two different forms: in Figures 1(b) and 1(c) the graph is symmetical under 
reflection in the vertical axis; this kind of symmetry is described by the formula 
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S(t) =f(—t). This is the kind of symmetry which also applies to the function 

cos wt, and a function with this symmetry is called an even function. On the other 
hand, in Figures 1(a) and 1(d) the graph is symmetrical under a rotation about the 
origin and this symmetry can be described by the formula f(t) = —f(—t) which 
also applies to the function sint. A function with this symmetry is called an odd 
function. 


Now it is a fact that the Fourier series for an even function does not contain any 
sine terms (i.e. A, = 0); apart from the constant term it can only consist of cosine 
terms. Similarly an odd function does not contain any cosine terms (i.e. B, = 0) 
but only sine terms. 


This piece of knowledge can save you the trouble of trying to evaluate coefficients 
which you can immediately predict will be zero. Look back at the four Fourier 
series in Section 2 and check that they conform with what I have said about odd 
and even functions. 


Exercise 1 


Figure 2 shows the graph of a periodic function of t in the shape of a triangular waveform. 
It is drawn so as to be symmetrical about the line ¢ = 0. It is an odd or an eyen function? 


[Solution on p. 25] 


If a function of ¢ has no symmetry about the origin then it is neither odd nor even 
and its Fourier series will contain both sine and cosine terms. 


Whether a function is odd or even or neither may depend merely on the choice of 


a St aa 9 S Figure 2, Triangular 
origin. For example, in Figure 2, if the graph were moved to the right by a waveform, 


distance less than a quarter of a period, the graph would no longer have symmetry 
about the origin and the Fourier series would contain both sine and cosine terms. 


Exercise 2 

State which of the quantities M, A,, B, will be zero for the function which results when, in 
Figure 2, the graph is moved 

(i) ‘one unit vertically downwards; 

(ii) ‘one unit vertically downwards and 4 period to the left. 

[Solution on p. 25] 


Exercise 3 
Derive the Fourier series for the function whose graph is shown in Figure 2, 


(You can take advantage of the symmetry when doing the integrations for the terms which 
are not zero by noting that, for example, 


2 2 
f S(t)dt = 2f f(thdt 
-2 lo 


because f(t) is even.) 
[Solution on p. 25) 


3.2 The use of Fourier series as approximations 


As we have seen, a Fourier series generally has an infinite number of terms. 
Clearly, in any numerical calculation we can cope with only a finite number of 
terms and the question therefore arises: how many terms do we need to get a 
good approximation to the original function? Of course, it all depends on what we 
mean by a ‘good’ approximation; and this must depend on what we want the 
series for. No general specification can be given, but some ‘feel’ can be obtained by 
looking at some functions like those we dealt with in Section 2 and comparing 
their graphs with those obtained by adding up an increasing number of terms of 
the Fourier series. You may be surprised at how close even quite a small number 
of terms will take us. 


By way of an example, we can start with the saw-tooth waveform, the first Fourier 
series we derived. Figure 3 shows the results of using from one to six terms of the 
series. 
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for KOA 
using one term using two terms 
- > 
t iL 
Koa OL. 
using three terms using four terms. 
- - 
t t 
KO fO* 
using five terms using six terms 
> - 
i t 


Figure 3. Graphs of partial sums of the terms of the Fourier series for the saw-tooth 
function. 


Figure 4 shows a triangular wave with, superimposed, the fundamental component 
(curve I) and the sum of the first two components (curve II) of the Fourier series. 
Neither curve is very far removed from the original function. 


fh 


Figure 4. Triangular wave with the first two partial sums of the Fourier series 
8 1 


ft) = (C08 + 530083 + 3 


a cos St + °°), 
Es 
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Clearly, then, even if a Fourier series is ‘truncated’ (i.e. cut off) after a finite 
ntimber of terms it can still give quite a reasonable approximation to the original 
function. 


Summary of Section 3 


1, A function which exhibits symmetry about the origin such that f(t) =f(—1) is 
known as an even function and its Fourier series contains no sine terms (A, = 0). 


2. If the symmetry about the origin is such that f(t) = —f(—t) the function is an 
odd function and its Fourier series contains no cosine terms (B, = 0), 


3. Fourier series are used as approximations to the original functions by using a 
limited number of terms, For many purposes good approximations can often be 
obtained with quite a modest number of terms. 


4 Extending the scope 


4.1 Functions of variables other than time 


So far the main applications we have had in mind have been those concerned with 
variations of time. It is quite possible, by analogy, to use Fourier series in cases 
where the independent variable represents a quantity other than time (e.g. space). 


Suppose we have a function of some other, unspecified variable, x. Suppose further 
that the function repeats itself at intervals of length L along the x-axis, as shown 
in Figure 1, We can derive the Fourier series exactly as we did before, provided 
that we substitute x for ¢ and L for T. We shall then finish up with a Fourier 
series in x, as follows, 


Fourier series in x 


f(x) =M+ ¥ (A,sinnaox + B,cosnwx) 
net 


Fe 2 2 
where is given by — = L, or » = o 
o 


1 pele 
and Mieos() ree 
L dein 


2 phe 

A, = = JS (x) sin nex dx 
L Jai 
2 phi 

B, ==> S (x) cos nex dx. 
L diye 


4.2 Some non-periodic functions 


So far, we have dealt only with periodic functions. It is possible, however, to 
extend the techniques of Fourier series to some non-periodic functions. This aspect 
of the subject is important mainly as a mathematical technique and this is how it 
will be used in the next unit. It is of very limited use for solving the kind of 
‘frequency response’ problem which we discussed in Section 1 because that kind of 
problem assumes a periodic input. The fact that I refer to ‘some’ non-periodic 
functions suggests that there is a restriction on the type of function we can cope 
with. The restriction is that the function must have a finite domain: it is defined 
over only a finite part of the t- or the x-axis. 


Figure 1, Periodic function of x. 
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The procedure for finding a Fourier series to represent such a function really 
amounts to extending the definition of the function to the whole of the f- or x-axis 
in such a way that the extended function is periodic. Then we can find the Fourier 
series for the extended function but use it only over the range for which the given 
function is defined. Let us take a simple example, the function F(x) shown in 
Figure 2. 


Here F(x) = 2 — x and the domain consists of the interval of the x-axis between 0 
and 2. The trick is to incorporate this function in a periodic function for which the 
Fourier series can easily be found. One way of doing this is shown in Figure 3. 


I have drawn dashed lines to define an even function (x) which is periodic with a 
period 4. Having got the Fourier series, we must remember to use it only over the 
range, 0 < x < 2, of the actual non-periodic function we started with. The Fourier 
series 1s like the one we derived in Exercise 3 of Subsection 3.1, Substituting x for 
t, it is 
2 8 Oe le aex. IL Sax 

Ti is Bite (C08 SOF == fC ) (for all x). 
This can be used to represent the function F(x) shown in Figure 2 provided that 
we restrict ourselves to the domain 0 < x < 2, ie. 


8 mx Le sax: 2 Sex 
F(x) =2 -— = 1+ sy (cos> + pc0s—- + se cos +++) (O<x <2). 
Exercise 1 


Using the first four terms of the series derived above, investigate how accurately this partial 
sum represents 


F(x) 


by substituting x=0, x=4, x=1, 
[Solution on p. 25) 


x forO0<x<2 


Exercise 2 


Again using only the first four terms, verify that the formula for F(x) given above cannot be 
used to represent the function 2 — x when x = 3, 


[Solution on p. 25} 


Exercise 3 


Sketch the graph of a periodic function which could be used to represent F(x) as defined in 
Figure 2 by means of a series containing only sine terms. (Do not calculate the coefficients.) 


[Solution on p. 25) 


Exercise 4 


Derive a Fourier cosine series with M = 4 and period 4 to represent the non-periodic 
function shown in Figure 4, for the range 0 <1¢ <1. 


[Solution on p. 26) 


Summary of Section 4 


1. The procedure for determining Fourier series for functions of time described in 
earlier sections can be used for functions of any variable x, with a period L, 
subject to substituting x for ¢ and L for T. 


2. A non-periodic function whose domain is a finite interval can be represented 
by a Fourier series by deriving the series for a periodic function of which the given 
function forms a part and then restricting the independent variable to the range of 
values of the variable for which the non-periodic function is defined. 


fit) 


Figure 4 
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5 End of unit problems 


Please refer to the study guide for this unit before you begin work on this section. 


Problem 1 
Derive the Fourier series for the periodic function one complete cycle of which is shown in 
Figure 1. 
Kos 
f}—2d—4 
hd ——— 


| 
1 
| 
1 
i 
_L 


Figure | 


[Solution on p. 26) 


Problem 2 
Figure 2 represents the frequency response of an idealized device, 


output amplitude 
input amplitude a 


- 
i 
1 w 
1 od 
| 
1 
i 
! 
1 
\ 
1 
i 
| 
Figure 2 
Show that, in response to any periodic input, the output of the device will be an exact fx), 
scaled version of the input, shifted by an amount t along the time axis. 
[Solution on p. 26} d 
Problem 3 


(i) Derive a Fourier sine series of period 2/ to represent the non-periodic function whose 
graph is shown in Figure 3 over the range 0 < x < |. 


(ii) Assuming that we may neglect terms in the series having an amplitude equal to or less 
than 1% of the amplitude of the first term, how many terms of the Fourier series are 
required to represent /(x) to this degree of approximation? 0: 


Solution on p. 27 
t a) Figure 3 
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Appendix 1: Solutions to the exercises 


Solutions to the exercises in Section 1 


1. (i) The fundamental angular frequency is 2 radians per 
second, 


(ii) The period is 2 = nseconds. 


2. (i) Since the angular frequencies of all the components 
are integral multiples of the lowest one, ic. w = 4, the period, 
T, of the function f(t) will be 2n/4 = n/2 seconds. 


(ii) The first component has angular frequency w = 4. At 
this frequency the relative response, or amplitude ratio, is 
found from Figure 7 to be about 1.2, and the phase shift 
about —11°, or —0.19 radians. Hence the output 
corresponding to the term sin 4t is about 1.2 sin (4r — 0.19), 


Similarly at @ = 8, the amplitude ratio is about 2 and the 
phase shift about —44° or —0.77 radians. The output 
component corresponding to the term 2sin 81 is about 

2 * 2sin (8¢ — 0.77) = 4sin (8¢ — 0.77), 


Again, at @ = 12, the amplitude ratio is about 1.5 and the 
phase shift is about — 132° or —2.3 radians, The component 
of the output corresponding to the term 3 sin 12r is about 
3x LSsin (12r — 2.3), Le, 4.5sin (12¢ — 2,3), 


Hence, by the superposition principle, the total response of 
the given linear system to the input f(t) is 


in (4t — 0,19) + 4sin (8t — 0.77) 
sin (12 — 2.3), 


x(t) 


de 


(Your answers may differ slightly because the graph is hard 
to read accurately.) 


Solutions to the exercises in Section 2 


n2 ra 
1, f J (t)cos cot dt = Mcosotdt 
2 


“Ta ae 


/2 
+ f A, sin wt cos ct dt 
=1/2 


T/2 
+f B, cos? wtdt + °°". 
-12 3 


All the terms on the right-hand side except that containing 
B, are equal to zero. Hence 


12 1/2 T 

f S(t)cosctdt = f B, cos? wtdt = By x — 

=n =T7/2 2 
Therefore 


2 pra 
By = =) f(t)cos corde. 
T deny. 


2. (i) The period would still be 2x, and so w = 1 as 
before. Also, from the symmetry of the graph, M = 0 as 
before. So the new Fourier series has the form 


S(t) = ¥ (A, sinnt + B, cos nt). 
m1 
Now, for —x <1 <2, f(t) = 2t/n, so 
PHYS Bihan 
A,= TF i i 7 ntdt = Fl if hi nt dt, 


This is just twice the original value'so 


Also 
Deze 
B= / = dt 
Bias tl Weta 


and again, since this is just twice the original value, B, = 0. 
So the new Fourier series is just twice the previous one: 
4/. sin2e | sin3r sin 4t 

fi-=4|[sine 58% 4 SB _ sind ) 
(ii) The period is doubled, so T = 4x and w = $. The 
Fourier series now has the form 

2 it 
S)=M+ ¥ (A,sin™ + Bycos”™), 
nt 2 2 

For —2n <t < 2m, f(t) =t/2n, so 


ed alraces 
a —dt =.—,| — =0 
Mee ee als. 


(as expected, since the mean value is still zero). 


1 pe t 2m 
A= x I xsinede= ae tsin dt 
LN (apo (aga 
@ rir ty 
= Gaig (~2xcosnn + 0 ~ 2c0s (nn) ~ 0} 
Sarin 
m 
1 pn 
B,= toos —dt 


So the new series is 
| HL eile } 


so=2 


r sin (t/2) - —- 


and the only difference is the changed value of w. 


3. From the graph, the period is 2. So m = 2n/T =n. 
Here, 


r/2 1 42 
(dt = t)dt = Ide 
ee a 


since over the rest of the cycle f(t) = 0. This time the mean 
value is clearly not equal to zero since no part of the wave- 
form is below the axis, So 


Lie Ca i 1 
== t)dt == ldt =| + — 
Pilewia am eeu 2\2 02 
7/2 1/2 
S(t)sin nwt dt = 1 
12 - 


il 


=> 


2 
Air 


ae 
ne J a/2 


1sin nat dt 
12 


" 
\ 
g|- 

EI 
— 
g 
1 
8 
oe 
eau 
aan; 
i 
© 


7/2 1/2 
{ F(0)e0smat de = [ cos nat dt 


a 

Tdera 1a 

-(3 sinm | = | sin —in =) 
nn ga rele 2 


Baines =0 


2n 
3n 2 
B=; as Bee 
2. 
By= 7 sin2dx=0 
B in d 
seein = and so on. 


Since the values of B, are the coefficients of the cosine terms, 
the required Fourier series is; 


1.2 cos3nt  cosSnt cos 7nt 
f@=5+= cosixt — ——— + —_— - —_— + 


The brackets contain an infinite series of cosine terms (the 
sine terms all being Zero), and again, the higher the 
frequency the lower the amplitude, 


4. From Figure 3, T= 4= 2 


ax 
Hence GT = = 
P 1 pte 
The mean_value = M = a) f(Ode 
T dona 
Wie py 
“aLf*"als], 


ape? 
A, = if Psinnar dt. 
4d 
The integral tables in the Handbook give 
2 
fesinarar = (- + FA cos at + 2 int AG, 
a @ a 
Here, a = nw so 


Sih ecrsaeees 42 
A=; = zia3) cosmor | : 


“2 
in net 
-2 


- ie + jag) 8(—2n0| 


1 5 e 
ae [2sin 2nw — (—2sin(—2nw))] 
=0. 


1pte 

Also B, = 3f ? cosnast dt. 
2), 

From the Handbook 


Je cosarat = 


BN 2t 
— 3] sinat + — cos at + C, 
a a 
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Now a = nw, as before, hence 
| en el) epee ia 
"31 Naw — na A 


BAe 


2 (+- 323) sin(—2n0) | 
+ ug [2e082n0 — (—2e0s(—2na))] 


But w = 7/2, so that sin 2nw = sinnnz = 0, and 


waga (I 
16 ‘e 
B= ax (“=a 
B,=- x (-1? => 
16 i ne 
Bs = 53 * (-1)= 55 and so on. 
So 
so=4 ug f= noose 
ok We taker re Chaar} 
4, 6s ate 
“ot 2 


5. From Figure 4 the period, T = 6h. 
So w = 2n/6h = x/3h. 


M=7] sod al Sloat 


ails ae hdt + ike tdt+ f “na) 


eo ay Be 
=a! 2h? + 0 + 2h?) =0. 


(In fact we can see directly from the graph that the mean 
value is zero.) 


7/2 3h 
A== a S(t) sin neot dt = af (0) sin neot dt 
7,2 3 


h o 
=all ery: tsin nest dt 
nh \ J 55 “hn 


rah 


+ i) hsin not a 


h 


Using the integral tables in the Handbook to find the second 
integral 


h 


1 t 1 
+ >| ——cos not + ——;sinnot 
al no war ] 


1 1 3h 
+ { —— cos nor 
3 nw 
—— [cos (—nwh) — cos (—3nwh) — cos 3nwh + cos nah] 
3ner 
1 


* Seok 


-" 


[- heosnwh + — sinneoh — heos (—nwh) 
no 


ber 
~ da 
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The terms in cos nwh cancel out, leaving 


Slt sinn@h — hcos ano. 


"3noh 
But @ = n/3h, so 
2 


~ hos) 


3h 
J (0) cos nest dt 
30 


-a(f = heosnandt + [ toos nent dt 
h\ Js, -h 


3h 


+} heosnot at). 
h 


Using the integral in the Handbook for the second integral, 


- 
B,= {-asiomor] 


1 
+a snot -+ soo ms | 


ae 3h 
Cy [sem ror 
3| nw hh 


= = [= sin(—ncoh) + sin (—3ncah) + sin 3neoh—sin noah) 
Ind 


oh 


: smal si noh + Leos no + hsin (—nwh) 
* Sncoh 


er 
no 


=0. 
Hence M'= B, = 0 and 
2h | 3 
== (= sin™ —cosnn 
nn \nt 


Solutions to the exercises in Section 3 


1. The function is even since f(—t) = f(t), and its Fourier 
series would consist of a constant term (since its mean value 
is non-zero) and cosine terms. 


2 (i) AO 


This will be an even function with zero mean value ie. 
M=4,=0,  B,#0. 


(ii) Ko 
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This will be an odd function with zero mean value ic. 
M=B,=0, 4,40. 
3, From the graph, T = 4, so w = 2n/4 = n/2. 
14 Bitton. 
ee 7) floa = Fala 
Between 0 and 2, f(t) =2 


1 (7 1 i 1 
M=5[@-od=3[ar-F] = ha -aa1 


Since the given function is even, A, = 0. 


1,80 


= 22 21008 dt 

a [ott dt 

2x2 sin SEE OL ald 
i om 2 ee 


4 4 1 ri 
= a. — sinnz — —cosnx — sind +0 + woos} 
nt nn nn 


4 
= wel —cosnm), 
8 
Thus B, = =o B,=0, By= om and so on. 
Hence the required Fourier series is 
SQ = 1+ s co! anes ot Pras 
ml anon wea ase o 


Solutions to the exercises in Section 4 


‘| 
1. F(0)=2 and 14S (145 +3 1.933. 


FQ) = th and 1+ S701 — 0.07857 — 0.02828) 


= 1.487. 
F(\)=1 and 145 0+0+0=1. 
ia 
S35, 
So that even with only a few terms the series gives an 


approximation to F(x) = 2 — x with a maximum absolute 
error of 0,067 for the stated values of x. 


F(2)=0 and 1+ 5-1 } = 0067 


8 
2,2-3=-1 and +=30+0+0)=1. 


This confirms that outside the stated domain the series does 
not represent the original function. 


3. The required graph is shown below. 
Se) 


Note that the full line represents the function G(x) = —x, 
0 < x <2, so that 2 would need to be added to the sine 
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series representing this graph in order to obtain the 
required representation of the function F(x) = 2 — x, 
O<x<2. 


4. We need to incorporate the function shown in Figure 4 
in a periodic function whose Fourier series will then 
represent the given function over its range of values of ¢, i.e. 
O<r<t. 


Exercise 3 in Section 2 provides a suitable function which, 
with a little adaptation, will give us what we want, We first 


fit) 


0 eae 
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draw the appropriate diagram in which only the given 
function appears in full lines. 


This is like Figure 2 in Section 2, except that the period is 
different. Here T = 4, so that w = 2/2 whereas in Section 2 
T =2 and w = x. This, however, is the only difference so 
that we can use the same series as we obtained then, 
changing only the frequencies of the terms. We therefore put 


AO heeee ost Loner 
Os (iat Tale aaa 


1) St Le 
+ Zcos-t — =cos—r + +++ 


3085 7085 forO<1<1 


You can check by substitution that this series works over the 
stated range 0 <t < |. With a ‘jump’ discontinuity at 1 = 1, 
the series will give the average value of /(t) for the vertical 
portion, ie. for ¢ = 1 the series gives f(r) = 4. 


Appendix 2: Solutions to the end of unit problems 


1, From the graph in Figure 1, T = 2I, so w = 2n/2l = n/l. 
The graph is symmetric about the t-axis so M = 0, It is an 
odd function so the Fourier series has sine terms only, i.e. 
B, = 0. 


2 p72 1 


f 
A= flo)sinnor de = + f f(t)sin a at 
T Jaa Wei U 

1 1 
a Nea 
=i(f = sin ae +f isin ™ at) 
Tee, 1 at 


= foto -2-4 
coho) ool] 
Amel A -ml 


Using the trigonometric identity from the Handbook 
sing sin B = 4cos (« — B) —4cos(x + f), 


— 4h ont. nnd _ nt 
S(= Die orn tar 
2. Let the input be f(t) where 
S(t) = M + A, sinat + Az sin Zot + Ay sin3at + °° 
+ B, coswt + B, cos 2wt + By cos 3a + + 


Now the amplitude of each of these terms will be multiplied 
by y in passing through the device since the same multiplier 
applies to all frequencies, (For this purpose M counts as a 
cosine term with zero frequency.) 


Each term will also experience a phase shift given by 
o= —to. 

For M, du =t x 0=0, 

for A, sinwt, @, = —ta, 


for A; sin er, o: = —2tm and so on, 


and similarly for the cosine terms. 
Hence, if the output is x(t), then 


x(t) = M+ WA, sin (wt — wt) + WA2 sin (Zot — 2wt) + +** 
+ WB, cos(wt — wt) + wB, cos (2mt — 2at) + -"* 
=p(M + A, sinw(t — t) + Az sin 2@(t — t) + -** 
+ By cosco(t — t) + Bz cos 2m(t — t) + **"), 
This is an exact scaled version of the input f(t) where (t — t) 
is substituted for ¢, ic, each term is shifted by an amount ¢ 
along the t-axis. The ideal frequency response shown in 
Figure 2 of Section 5 is not, in practice, attainable. You may 


be interested to see how near to the ideal a practical 
measuring device in common use gets. 


amplitude ratio (%) 4 


0 rf 
0 02 04 0.6] 08 10 1.2 14 


frequency ratio = applied frequency 
natural frequency 


frequency ratio 
02 04 06 O08 10 12 14 


(radians) 4 
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The above figure shows the frequency response graphs for a 
galvanometer on a commercially available voltage recording 
instrument. You can see that up to a frequency of about 0.7 
times the natural frequency of the galvanometer, the actual 
relative response is within +5°% of the ideal, and the phase 
shift graph is quite close, too, For frequencies higher than 
about (0.7 x natural frequency) the relative response departs 
very markedly from the ideal. Hence, in order to get a true 
recording of a periodic input from this instrument, the terms 
of the Fourier series of the input with frequencies higher than 
(0.7 x natural frequency of galvanometer) must be negligible 
or zero. 


3. (i) The question asks for a sine series, ic. a series 
consisting only of sine terms; this means that we must find 
an odd function of which the given figure is a part. Such a 
function is shown below, 


feo 


ay 


\/ 
be | 


Here the ‘fictitious’ parts of the periodic function are drawn 
in dotted lines. For this periodic function, M = B, = 0, and 
L=121, so that w = 2n/L = x/l. 
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But » =" 50 sin nol = 0 and 


T 
8d on 
A.= 780 
8d 8d 8d 
Thus Ay =, 42 =0, dy = ~ 5 Ae = 0, As = 555 and 


so on. 


Hence the Fourier sine series is 


Re (en an gi 
ie) eT a 
for0<x<! 
or 
8d & (-1)"! n— 1ymx 
L(x) Beit for0<x <I. 


(ii) The amplitude of the first term is “ 


The ratio of the amplitude of the nth term to that of the first 


term is @-1 


. 1 ie 
For this to be greater than 1% we want @n-1F > 100 ie. 
(2n — 1)? < 100 or 2n — 1 < 10. This means for n < 5 we 
must include the terms in the Fourier series, so five terms are 
required. 
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Appendix 3: Television programme notes 


You should read these notes before watching the television programme ‘Fourier 
Analysis’. 


A.1 Introduction 


The television programme for this unit is about measurement and recording. It 
describes a procedure for predicting the response of a recording device to any 
periodic input that can be represented by a Fourier series. This procedure can be 
used for instruments such as audio amplifiers, record players or 
electrocardiographs. In the programme we use it to test a galvanometer which 
forms part of an ultra violet light recorder. 


In order to use the procedure we need two pieces of information: 
(i) the frequency response of the recording device; 
(ii) the Fourier series of the periodic signal to be measured or recorded. 


Since the Fourier series has an infinite number of terms it must be truncated to be 
used in practice, but enough terms must be left to constitute a good 
approximation to the original waveform. As I said in Subsection 3.2, there is no 
hard and fast definition of a ‘good’ approximation; the programme models the 
triangular wave, which serves as the main example, as the sum of only the first 
three terms of the relevant Fourier series. (You can get an idea of how good this 
approximation is by looking at Figure 4 in Subsection 3.2 which shows the sum of 
the first rwo terms in this Fourier series compared with the complete waveform.) 


The programme also shows how the relative response of the galvanometer is found, 
by providing a constant-amplitude sinusoidal input over a range of frequencies and 
noting the variation in the output amplitude. 


A.2 A note on frequency 


In the programme, frequencies are quoted in cycles per second because 
instruments are calibrated in this way. However, in this unit we have worked 
exclusively in angular frequencies (radians per second). The connection between 
these quantities is quite simple. 


If J = frequency, in cycles per second (or Hz) 
and  =angular frequency, in radians per second, 


then 
and 

w = 2nf. 
Example 1 


(i) The fundamental frequency of a periodic waveform is 30 cycles per second. 
Calculate the corresponding value of the angular frequency. 


(ii) The component with highest frequency in a truncated Fourier series has an 
angular frequency of 600 rads~'. What is its frequency in Hz? 
Solution 
(i) @ =2nf~ 188.5 (in radians per second). 
@ 600 


(ii) f= of ae ~9549 (in Hz). 


The relative response, i.e. 
output amplitude 
input amplitude * 

is called the amplitude 

response in the television 

programme. 


The SI unit of frequency is the 
hertz: 
1Hz = 1 cycles™!. 


fand o are related to the 

period, T (in seconds), by 
ple eS 
ay Oma 
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A.3_ Relative response and phase shift 


As you saw in Subsection 1.2, the frequency response of an instrument consists of 
its relative response (i.e. output amplitude/input amplitude) and its phase shift. 


The relative response of the galvanometer in the ultra violet light recorder is 
shown in Figure 1(a). For good reproduction, the graph of the relative response 
should be as nearly as possible parallel to the horizontal axis over a range of 
frequencies. This range should include all the frequencies in the truncated Fourier 
series representing the signal to be recorded. In this case you can see that the 
graph is flat up to a frequency of about 90 cycles per second. 


output amplitude 4 


input amplitude 


100” frequency 
(a) {eyeles per second) 


0 


frequency 


phase shift 4 
(cycles per second) 


0 


(b) 


Figure 1, The frequency response of the galvanometer tested in the television programme. 


The phase shift arises because the recording device takes time to respond to an 
input signal so that the output is not simultaneous with the input; each 
component of the Fourier series representing the input is delayed by an amount 
which depends on its frequency. Even if all the output amplitudes are correct, they 
will not add up to a good reproduction of the input if the instrument introduces a 
time shift of one component relative to another. In other words, the time shift 
should be the same for all components. To see what this means in terms of the 
frequency response, consider just one sinusoidal component, f; = F;sin wt, say, of 
the input, f(t). Then the corresponding component of the output, x(t), will be 

x, = X,sin w(t — t), where t is the delay introduced by the instrument. But in 
terms of the phase shift, @, the output component is x; = X,sin (wt + ) so 

= -or. (1) 
We want the time shift, z, to be constant for all values of w, so Equation (1) shows 
that the graph of @ against w should be a straight line. The relevant graph for the 
galvanometer is shown in Figure 1(b), and you can see that the frequency range 
over which this graph approximates a straight line is greater than the frequency 
range for which the relative response graph is flat. So provided we stay within the 
frequency limits imposed by the relative response, we shall automatically satisfy 
the phase shift requirements as well. For this reason the programme concentrates 
on the relative response. 
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Example 2 

Assume that the triangular waveform shown in Figure 2 can be represented by the 
sum of the first three terms of its Fourier series; estimate the highest fundamental 
frequency at which such a wave will be reproduced by the galvanometer used in 
the programme. 

Solution 

From the solution of Exercise 2(i) in Section 3, M = A, = 0, Modifying the 
solution of Exercise 3 in Section 3, we find that the series is 


8K 1 1 
So®= gz (cost + 9008 3a + 75 008 Set) 
where «w = 2n/T is the fundamental angular frequency and K is the height of the 
triangular wave. 


Now, in order to stay within the flat region of the relative response graph (and 
~ hence also the linear region of the phase shift graph), we must have 


5m < 2m x 90 


2n x 90 
5 
The fundamental frequency must be less than about 113 radians per second. 


as 


e118; 


Now watch the television programme ‘Fourier Analysis’. 


